Arrays of Josephson junctions between unconventional superconductors 
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We study large arrays of mesoscopic junctions between gapless superconductors where the tunnel- 
ing processes of both, particle-hole and Cooper, pairs give rise to a strongly retarded effective action 
which, contrary to the standard case, can not be readily characterized in terms of a local Josephson 
energy. This complexity is expected to arise in, e.g., the grain boundary and c-axis junctions in 
layered high- Tc superconductors. A new representation for describing collective phenomena in this 
system is introduced, and its phase diagram is discussed, alongside the electrical conductivity. 



Quantum dynamics of ultrasmall normal and super- 
conducting (Josephson) junctions (JJ) has long been a 
field of active theoretical^ and experimental^ research. 
Recently, the interest in this topic has been rekindled 
by a number of new ideas, such as the proposal of a 
novel 'floating' phase, in which context the effects of 
(spatially) long-range correlations were investigated at 
a greater length^. 

Notably, though, most of the previous theoretical stud- 
ies were limited to the JJs between conventional, fully 
gapped, s-wave superconductors. Although the case of 
the d-wave superconducting cuprates, such as bi-epitaxial 
grain boundary (in-plain) JJs in YBCO or intrinsic c- 
axis (vertical) ones in Bi2212, have been rather exten- 
sively studied as well, their previous analyses would rou- 



tinely resort to a phenomenological description similar to 
that of the gapped (s-wave) superconductors^. In con- 
trast, the microscopic analysis of a single d-wave J J car- 
ried out in Refs.- t'', showed that the processes of both, 
particle-hole and Cooper, pair tunneling can give rise to 
the equally non-local (in the time domain) terms in the 
effective action, thereby invalidating the very notion of a 
local Josephson energy. 

In the present work, we study a large array of such JJs 
and look into the effects of a strong time dependence of 
the effective action on this system's phase diagram and 
corresponding behaviors. 

The partition function of a generic single JJ introduced 
in Ref.^ can be generalized to the case of an array by 
including both, self- and mutual-, capacitances 
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where (r) = (pi {T)~'(j)j (r) is the phase difference across 
the link < ij >. 

The double time integrals in (1) are governed by the 
kernels a(T) and (3{t) representing particle-hole and 
Cooper pair tunneling processes, respectively. To the 
leading order in the tunneling matrix element T{k,k'), 
they are given by the expressions 
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where Q and T are the normal and anomalous electron 
Green functions, respectively. 

The a-term describes (non-Gaussian) dissipation due 
to the Andreev quasiparticle tunneling whose effects 
have been extensively discussed in the previous worksi, 
while the /3-term represents the processes of (in general, 
non-synchronous) pair tunneling. In the conventional 
(s-wave) superconductors, it decays as /3(t) oc e^'^l'^l, 
thereby effectively reducing the last term in (I) to a 



single time integral Ej J^^^ dr cos20ij(T) of what can 
then be identified as the local Josephson energy E,j = 

By contrast, in the case of a gapless superconductor 
one obtains strongly retarded kernels^ 

a{r)/a - /?(r)//3 - l/r^^"" (3) 

where the prefactor in the /3-kernel vanishes for any fac- 
torizable matrix element, T{k, k') — f{k)f{k'), of a sym- 
metry other than the s — wave and, therefore, it can only 
be due to a non-factorizable contribution oc f{k~k') into 
T{k,k'). 

In the two-dimensional case and under the condition 
of momentum conservation, 2D — rj = 2 (see Refs.^'^), 
hence both the tunneling terms appear to be marginal, 
the corresponding coupling constants a and /3 being di- 
mensionless numbers of order unity. A short-time diver- 
gence of Eqs.(3) can be naturally regularized by substi- 
tuting T Vt^ + A~2 where the cutoff scale A is set by 
the maximal superconducting gap in the bulk. 
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Conceivably, one can encounter even longer-ranged 
correlations (2D— 77 < 2) in the presence of, e.g., resonant 
tunneling through zero energy states supported by cer- 
tain tunneling configurations, such as that of the (io/'^Tr/4 
in-plane grain boundary^. 

Turning now to the effective action (1), we find that 
a strongly retarded nature of the tunneling terms ren- 
ders a customary dual representation based on the Vil- 
lain transformation of the local Josephson term inappli- 
cable, thereby making this model unsuitable for the stan- 
dard mapping onto an effective vortex plasmEti. There- 
fore, a well-known description of the different phases 
in terms of bound vortex-antivortex complexes (dipoles, 
quadrupoles, etc.) can not be readily generalized to the 



problem at hand, either, thus forcing one to take a dif- 
ferent approach. 

To that end, we introduce a new bosonic field ipiir), 
alongside an associated Lagrange multiplier field enforc- 
ing the local constraint i^iir) — e"^'^'^-'. This approach 
should be contrasted with the previously developed treat- 
ments of the conventional (local) Josephson term (see, 
e.g., Ref.*) where a constrained bosonic variable would 
be used to represent the pair field e^"^'*^"^^. Indeed, an 
attempt to implement this technique in the present (non- 
local) case would require one to work with a technically 
intractable bi- local composite operator 'ipi{T)Tpi{T'). 

By integrating out the phase variable 0i, keeping the 
leading terms of the corresponding cluster expansion (cf. 
with Ref.^), and then integrating out the Lagrange mul- 
tiplier field, one arrives at the partition function 
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where Ai(r) is an additional Lagrange multiplier enforc- poratcd into the normal Gij =< V'iV'j- > ^^id anoma- 



ing the auxiliary constraint 'ipl{T)Tpi{T) = 1 (the latter 
is not automatically satisfied, unless the integration over 
(j>i{T) is performed exactly). 

The correlation function appearing in Eq.(4) 
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is governed by the effective Coulomb energy Ec = 
I 2(27i-'jp+iCfc proportional to the integral of the inverse 

capacitance Ck — J2<ij> Cije^''^^^ which converges, pro- 
vided that the capacitance matrix progressively decreases 
with the separation between the sites. 

The frequency integral in Eq.(5) diverges for any Rij ^ 
which dictates that the correlation function Wij (r) re- 
mains strictly local in the real space. Also, Eq.(5) is 
written in the limit of vanishing temperature, while at fi- 
nite T a proper account of large phase fluctuations with 
non-trivial winding numbers makes this (as well as any 
bosonic) function periodic with a period 1/T by virtue 
of the substitution t — > r — Tt^ (see Ref;i). 

At a = /? = one then obtains a bare (normal) Green 
function 
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while for finite a and /3 the quantum charge fluctua- 
tions give rise to the corrections which can be incor- 



lous Fij —< ipiipj > Green functions obeying the usual 
Dyson's equations 
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where both the normal and anomalous Ay self- 
energies can be computed as series expansions in powers 
of a and /3. 

The analysis of these expansions shows that they can 
be organized according to the powers of the inverse co- 
ordination number z (e.g., z = 2D for a simple cubic 
lattice). In the leading approximation for z ^ 1, the 
self-energies are given by the equations 
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When ascertaining a general layout of the phase diagram 
of the JJ array, different components of the self-energy 
can serve as emergent order parameters. As such, one can 
distinguish between the local, Eq = ^mj ^-nd non-local, 
Si = i (here the sum is taken over the z near- 

est neighbors), normal, as well the corresponding anoma- 
lous, Ao = Aii and Ai = i Aj,i+^, self-energies. 



Specifically, Si signals the onset of a metallic behav- 
ior (hopping between neighboring sites) , Ap manifests an 
incipient local pairing, Ai serves as the precursor of su- 
perconducting coherence setting in across the entire JJ 
network, while a frequency-dependent part of the Sq in- 
dicates the development of local time correlations. 

With the on-site and nearest-neighbor terms taken into 
account, the spatial Fourier harmonics read 
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where 7(fc) — t 

Eqs.(8) can be further improved by adding polarization 
corrections to the effective coupling terms 
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where the polarization functions IIe.o{^) = 
J ^Te,oG{uj')G{uj - uj') include the vertex cor- 
rections Te,o arising from the even and odd numbers of 
non-crossing /3-couplings 
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With the vertex and polarization corrections included 
and in the absence of any emergent order parameters, 
the self-consistent equation for Eo(a;) reads 



Lnlo; — Z i (OJ , 
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The (static and spatially uniform) expectation value 
of the Lagrange multiplier A = (Ai(T)) can then 
be determined from the normalization condition 

In order to ascertain the locations of the putative phase 
boundaries we include a constant term I]o(0) + A into 
the definition of the renormalized Coulomb energy Ec 
and expand Eqs.(8) to the first order in the emergent 
self-energies Ei, Aq, Ai, as well as the derivative of the 
(linear) frequency-dependent part of So(w). Threshold 
values of the couplings, beyond which such self-energy 
components develop, are then given by the eigenvalue 
equations 
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FIG. 1. Left panel: The onset of the inter-site self energy 
El and both on-site and inter-site anomalous self energies Aq 
and Ai. Right panel: phase diagram (see text) 



In the case of marginal ('Ohmic') dissipation correspond- 
ing to 2D — 77 = 2 the Fourier transforms of the (regular- 
ized) coupling functions behave as a{u!)/a = /3{uj)/l3 — 
TrAe^l'^l/'^, thus resulting in only a weak frequency de- 
pendence of the self-energy at w <C A. 

The first three of the eigenvalue equations (13) then 
reduce to the algebraic ones 
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from which one determines the locations of the putative 
critical lines in the a — /3 plane (see Figl). 

Interestingly enough, Eqs.(14) suggest that for small 
a and large z the onset of local ('on-site') pairing upon 
increasing /3 may precede that of the metallic behavior, 
while for small (3 the inter-site ('bond') pairing emerges 
only at sufficiently large a. 

The above observations suggest a general layout of the 
phase diagram presented in Fig.l. The region of small 
a and /? with Si = Aq — Ai = is interpreted as uni- 
formly insulating (I) , while the emergent order parameter 
Aq 7^ signals the onset of local superconducting pair- 
ing (LP) at /? ~ 1/z. At still higher values of /? ^ 1 
one expects to enter a Josephson-like phase (J) with 
Aq, Si 7^ but without global coherence. On the other 
hand, at a ^ 1 the insulator gives way to the resistive 
phase (R) with Si, Ai 7^ which supports both. Cooper 
pair and single quasiparticle, transport. Lastly, the uni- 
formly superconducting phase (SC) with Si,Ao.i 7^ 
would eventually be attained at a,/? 1. It should be 
noted, though, that our predictions are based on the ap- 
proximate perturbative analysis and, therefore, not all 
the putative phase boundaries may actually be present 
in the real system. In particular, there may or may not 
be a physical distinction other than a crossover between 
the J and LP phases, or the latter regime might be absent 
altogether (as it is for z = 2). 
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Such caveats notwithstanding, the overall behavior ap- 
pears to be somewhat reminiscent of that in the standard 
(s-wave) case: the system can be nudged closer to the su- 
perconducting state by increasing either, the Cooper pair 
or particle-hole tunneling, the latter providing a mecha- 
nism for intrinsic dissipation which quenches phase fluc- 
tuations and promotes the classical Josephson effect. 

Should, however, the tunneling /3-term happen to de- 
cay even more slowly {2D — rj < 1), the analog of the 

I 



effective Josephson energy would then diverge at large 
T, thus making the infrared behavior essentially singu- 
lar and possibly allowing for some drastic changes in the 
phase structure. 

Conducting properties of the JJ array allow one to 
discriminate between the different phases. In particu- 
lar, electrical conductivity can be computed as a^ii{ui) = 

'LtaT^ with the use of the action of Eq.(l) in the pres- 
ence of an external vector potential A^, resulting in 
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where the dots stand for 'paramagnetic' terms contain- 
ing higher powers of a and /3 which, therefore, are small 
compared to the above ('diamagnetic') contributions for 
a, /? ^ 1 (cf. with the discussion of a normal granular 
metal where /3 = in Ref.^). 

The thus-obtained longitudinal conductivity reads 
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and, upon performing the frequency integrations, one ob- 
tains 
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where, for the sake of simplicity, we chose T ^ Ec — A. 

The emergent metallicity order parameter Ei promotes 
a metal-like (temperature-independent at T — >■ 0) con- 
ductivity, thereby distinguishing it from the activation- 
type behavior characteristic of the insulating regime. In- 
terestingly enough, it also contributes to the superfluid 
density, alongside the local pairing Ag, while the non- 
local one (Ai) does not (to the lowest order in /?). 

It is conceivable, though, that there might be a (par- 
tial) cancellation between the 'diamagnetic' and 'para- 



I 

magnetic' terms at a, [3 ~ 1, as a result of which the 
conductivity could remain universal along the critical 
lines, akin to the situation in the conventional, s-wave, 
JJ networks^*' (it is worth reiterating that in the present 
case one can not readily invoke the charge- vortex duality 
on which the universality argument is based^ due to the 
inapplicability of the underlying Villain transformation) . 

To summarize, in the present work we studied arrays of 
unconventional JJs with long-range (in the time domain) 
interactions stemming from the presence of gapless quasi- 
particle excitations. On the technical side, the problem 
presents a new challenge by not being amenable to the 
customary approaches exploiting the intrinsic locality of 
the standard Josephson effective action. 

By using an alternative representation, we find that the 
phase diagram of the system might feature the insulat- 
ing, uniformly superconducting, Josephson (local pair- 
ing only), and metallic phases which can be identified 
by the corresponding emergent order parameters. We 
also predict that this picture might be further altered in 
the presence of resonant tunneling between zero energy 
states where the temporal decay of correlations can be 
even longer-ranged. 

We conclude by expressing a hope that this analysis 
will prompt a further investigation into (and provide an 
alternative means for interpreting the experimental data 
on) the assemblies of high-Tc JJs beyond the scope of the 
customary phenomenological approach adapted from the 
earlier studies of the s-wave superconductors. 
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